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Abstract 

We consider geodesies motion in a particular Kundt type III spacetime in which Einstein- Yang- 
Mills equations admit solutions. On a particular surface as constraint we project the geodesies into 
the (x, y) plane and treat the problem as a 2-dimensional one. Our numerical study shows that 
chaotic behavior emerges under reasonable conditions. 
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Kundt's class of solutions present spacetimes with non-expanding, shear-free and twist- 
free null geodesic congruences [l|. Interest in this class of solutions comes from the fact that 
these spacetimes admit plane waves, which exhibit geometrically different properties than 
the pp-w&ves. Due to itspotential application in string theory Kundt class still maintains 
its popularity Various Petrov types of Kundt solutions were identified, among 

them especially Kundt type N has been studied in detail for long time by many authors, (see 



for instance 3, 



IB 



). Geodesies in a specific Kundt type N was analyzed in detail by 



* , and i« was shown that particnlar so.ntions ohey a power-law. Recently, interest in 
the Kundt type III spaces has gained momentum. Firstly, Griffiths et al ]10[ derived and 
classified a complete family of Kundt type III, which admit cosmological constant A c and/or 
pure radiation $22 ^ 0. In the absence of the cosmological constant A c , the Kundt type 
III solutions are further generalized to be the solutions of the Einstein- Yang-Mills (EYM) 



system by Fuster and Holten [11]. 



The pioneering study of proving chaos in the spacetimes of plane waves was done for 
pure impulsive gravitational j?p-waves [.12] . Recently, it has been further demonstrated that 
under certain conditions the emergence of chaotic motion is possible both in the spacetimes 
of the superposed electrovac pp-waves and in the non-Abelian plane waves of Kundt type 
N [3], which are the solutions of the D = 4 EYM equations. In [ljj, it was also pointed 
out that for the non-Abelian plane waves the chaotic effect of gravity dominates over the 
gauge field. This is due to the fact that such local fields vanish asymptotically and the chaos 
inherited from gravity renders the whole system chaotic. 



Analysis of the geodesic motion for the D = 4 Kundt type III is also discussed by ll| in 
which a possible chaotic motion is highlighted in particular cases of this class. 

In this Brief Report, our aim is to investigate the possibility of chaotic geodesies in 



Kundt class II 
line-element 



spacetimes. For this purpose, we consider the following algebraically special 

u, 

ds 2 = 2du[dv + Hdu + Wdz + Wdz] - 2dzdz, (1) 



where H = H(u, v, z, ~z) is a real function while W = W(u, v, z, ~z) is a complex function, 
in general. Here our motivation is to study the geodesic motion for the case W >v = j^=, 
referring to the solution of the EYM equations for the metric (1) given by [ll| 
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2v 

W = W°(u, z) 



z + z' 



H-H'^ + ^l?™ ,-* (2) 

z + z (z + z) 2 

where W°(u,z) is an arbitrary complex function and H°(u,z,~z) is a real function. The 
simplest choice of W°(u, z) for the type III is 

W°(u,z) = g(u)z, (3) 

such that the Weyl scalar \l/ 3 7^ 0. On the other hand, imposing the only solution (x a = 
X a (u)z) on the Yang-Mills (YM) equation in which the energy density is bounded throughout 
the spacetime, the solution for H°(u,z,z) becomes 



H°(u,z,z) = [f(u,z) + f{u,z)] (z + z)-ggzz + cr(u)(z + z) 2 [Hz + z)-l]. (4) 

where f(u, z) is an arbitrary complex function and a(u) = 2 / y ab X a (u)\ b (u) + g(u)g(u) is a 
real function. Here, X a (u) are bounded complex functions and 7^ is the invariant metric of 
the Lie group. We note that the condition for being Kundt type III spacetimes is g{u) ^ 0. 
It is trivially seen that X a (u) = corresponds to the vacuum solution. 

Our primary interest here is to write the geodesies equations for the metric (1). Similar 
to the study III], for the beginning, we eliminate the w-dependence from W° and g by the 
following particular choice 



and g = 1, 



(5) 



Next, introducing real spatial coordinates x and y by \/2z — (x + iy), we get the geodesic 
equations as 



2 v x 
u-u (1 - — ) + 2u- = 0, 



X 



(6) 



. 2 

X + u 



H, r - 1- 



x 



X 



XJ J 



2 . . . . / v 
+ -uv + 2ux 1 - 2 — 



x 



(7) 



y + u 



H y + y ( 



X 2 



2ux— 

x 



(8) 
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where the dot denotes -f- with r being the proper time. In addition the metric condition 



implies 



.2.2 .. .2 .. V 

x + y — 2uv — 2Hu — 2ux{x — 2— ) + 2yuy = e. (9) 



where e = 1,0,-1 stands for timelike, null and spacelike geodesies, respectively. The 
present form of the equation set does not allow us to obtain a 2D (x, y) Hamiltonian system 
analogous to the previous studies [3, Q|- However, with appropriate choices of u and v 
surfaces it is possible to project the geodesies into the (x, y) plane in which writing a 2D 
Hamiltonian becomes possible. Our first intention is to shift the independent variable from 
r to u as an affine parameter. If we consider a family of geodesies, which follow the light-cone 
coordinate u with constant rate of change in the same proper time intervals, the following 
assumption can be made 

u = constant = 1, (10) 
Such an assumption gives rise to a condition on the v surfaces given by 

v = x\l~2-), (11) 

Here " ' " denotes By this substitution into Eqs. (7) and (8), we get a 2D dynamical 
system in the (x, y) plane 

3x" - H° x - x = 0, (12) 

y" + H° y = 0, (13) 
which is described by a Super-Hamiltonian [14] 

(14) 



2 V y 3 

with the corresponding potential 



V(x,y) = H° + j. (15) 



4 



Let us note that the Super-Hamiltonian defined by the momenta P x = — 3x' and P y = y' 
is not positive definite. 

Eq. (9) stands for an energy condition in which it should be automatically satisfied by 
the solutions of Eqs. (12) and (13). Without loss of generality, we can assume that / and 
A a are independent of u. This assumption implies that a is a positive constant. As we 
mentioned before that the chaotic effect of gravity dominates over the gauge (YM) field [3] , 
it would be sufficient to investigate chaos in vacuum, i.e. a = 1. In other words, once the 
chaotic motion appears in the vacuum spacetime, the local fields could not be strong enough 
to avert it into a regular motion. 

In general, any / = kz n (n = 0, 1, 2...), with the multiplicative factor k being an arbitrary 
real parameter implies a potential 



V = V2kxRe(z n ) - V — + 2x 2 \n(V2x) - 1 



(16) 



which admits an integrable system for k = 0. The logarithmic term in the potential 
imposes a condition on the x coordinate, namely x > 0. Beside this, the case / = kz, which 
describes a flat space for vacuum pp-wave spacetimes , confesses a regular motion for 

the geodesies particles. Contrary to the vacuum homogeneous pp-waves |l2j , here the case 
n = 2 admits a nonintegrable dynamical system. However, the nonintegrable systems in the 
vacuum pp-waves were obtained for the cases with n > 3, 12| • 

Here we wish to study the nonintegrable system with the simplest case (n = 2), and 
explore whether the motion depends on the initial conditions or not. If chaos emerges in 
such a simplest case, intuitively it should also appear for the / functions with higher powers 
of z. Now, for n — 2, it can be seen that the potential (16) has various unstable points 
according to the range of the multiplicative constant k. 

TABLE I: Unstable points depending on the value of k. 



Points 


Saddle 


Repellor 


V-i^V^ + SM 2fc) 


} 


k < -|e~£ 




jx = ^LambertW (ffc^e) , y = oj 




k ^ -§e~s 


4 a 7 1 7 

— |e 2 < k < — s 



Those results in the table show us that the highest possibility of the emergence of chaos 
corresponds to the case k < — \e~i , in which admits two saddle points. Particularly, 



the case 



\e 2 



< k < 



-,e s 



has an additive repellor point, and whence it may follow a 
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stronger chaos. Conversely, the case k ^ — \e~i upon possessing one saddle point causes a 
questionable chaotic motion. In order to judge the existence of the chaotic motion, we study 
the numerical analysis of the evolution of the test particle in the gravitational field. 

We integrate numerically the equations of motion given by Eqs. (12) and (13). The 
initial conditions depend on 3 parameters, (x , yo) (at u — 0) and k. For a given k value, 
we may choose (xq, yo) such that keeping xq unchanged, and checking the effect of the yo on 
the geodesic motions, while it varies. To do this, we may set x — Ci, a real constant c\ > 0, 
and yo — — 3+X/|- For example, if we take k = 0, the solutions are trivially analytic. This 

j=0 

is also graphically verified in Fig. 1. Next, by considering the cases k ^ the motion can 
lead to a chaotic motion. It is observed that the chaos has a movable character depending 
on the choices of Xq and yo while k ^ — |e~i, see Fig. 2. However, when the multiplicative 
constant k < — |e~i chaos is obvious. It is seen that the multiplicative constant k of the 
function / becomes decisive for the chaotic motion. In other words, k plays the role of 
critical parameter for the onset of chaos. The chaotic behavior of the geodesies in the case 
k < — |e~s is illustrated in Fig. 3. Alternatively, the dynamical system can be investigated 
by using the Poincare section method. We use the package POINCARE [14] to perform the 
numerical experiments. Fig. (4) is a demonstration of the Poincare section, which verifies 
the chaotic behavior in our dynamical system. 

In conclusion, it is shown that the Kundt type III spacetimes may reveal chaotic motion 
under certain conditions. To our knowledge such a study did not exist in the literature 
before. Chaos in the spacetimes of electrovac and the specific Kundt type N with YM field 



was studied before 



13J. This report constitutes an extension of that study. It is needless 



to state that the existence of the chaos in the Kundt type III spacetimes may have further 
implications for the particle motions in string theory and in higher dimensions. 
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Figures 



FIG. 1: For k = 0, 2D (x,y) plot of the geodesies. Geodesies start from 

18 } 

x = 0.7, i/o = — 3 + r (dashed line) and move through non-intersected trajectories. 

3=0 J 

The non-intersected trajectories represent the regular motion. 

FIG. 2: In case k > — |e~l, k — 1 is chosen for the 2D (x,y) plot of the geodesies. The 

f 18 1 

initial positions are < x = 0.7, ?/ — ~~ 3 + > (dashed line). Intersected trajectories sig- 

I i=° J 

nal the existence of chaos. Two trajectories, which have different yo initial points initially 

accelerate in +x-direction contrary to the others. 

FIG. 3: In case k < — |e~i, k = -0.21 (-|e"5 < -0.21 < -\e~l) is chosen 2D (x,y) 

r is.i 

plot of the geodesies. The initial positions are < x = 0.7, yo — — 3+ J^f } (dashed line). 

The chaotic behavior is evident from the trajectories. The symmetry in Eq. (13) shows 
itself along the y-axis. 

FIG. 4: Poincare sections of (x', x) for k = —0.21 and H = 0.2 across y = KAM surface. 
Some points are distributed randomly in a finite region to form a chaotic sea, however the 
large island surrounded by the chaotic sea indicates the existence of quasi-periodic orbits. 
(Here, x — > ql, and x' — > pi) 
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